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ABSTRACT 



Aims. We assess the possibility that planets Gliese 581 c and d are within the habitable zone. 

Methods. In analogy with our solar system, we use an empirical definition of the habitable zone. We include assumptions such as 
planetary climates, and atmospheric circulation on gravitationally locked synchronous rotation. 

Results. Based on the different scenarios, we argue that both planets in Gliese 581 could develop conditions for a habitable zone. 
Conclusions. In an Earth-like environment planet d could be within a habitable zone, if an atmosphere producing greenhouse effect 
of 100K could have developed. If the planets are gravitationally locked-in, planet c could develop atmospheric circulation that would 
allow it to reach temperatures consistent with the existence of surface liquid water, which in turn could support life. 

Key words, stars: individual Gliese 581 - stars: planetary systems super-Earth planets astrobiology 
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1. Introduction 

M stars are the most common stars in our galactic vicinity and 
about 75% (Tarter et al. 2006) of all stars in the Galaxy. Because 
of a lower mass (compared to our Sun) the M stars radiative 
output, after an initial highly variable stormy period, becomes 
stable for many billions of years, raising a possibility of evolu- 
tion of life on some of their planets. Following the discovery of a 
Neptune-mass planet around Gliese 581 (G58 1 hereafter), planet 
b, by Bonfils et al. (2005), Udry et al. (2007) discovered two ad- 
ditional super-Earth planets, G581 c and d, and, based on solely 
orbital considerations, determined that planet c was within the 
habitable zone (HZ). 

Later, adding atmospheric considerations or greenhouse ef- 
fects, von Bloh et al. (2007a) recently concluded that planet c is 
outside the HZ, whereas planet d is well within the HZ. 

G581 has a mass of about 0.3 IMq, spectral type M3 dwarf, 
and its photospheric temperature is estimated to be T e ff = 
3480/T, compared to the T eff of 5110K for the Sun. G58 1 lumi- 
nosity has been estimated to be Lc58i = 0.013 + 0.002Lq- The 
planets' masses and their distances (semi-major axis) from the 
central star are M c = 5.06M ffi , R c = 0.073 AU, and M d = 8.3M ffi 
with R f / = 0.25 AU, respectively. 

In this Letter, we consider the physical conditions of plan- 
ets c and d to be similar to conditions of the inner planets in 
the solar system and we show that under special, but feasible 
circumstances, both planets could have developed HZs. 
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2. Limits on the solar system HZ 

HZ is generally defined as a zone within which a planet can 
maintain liquid water on its surface. Different sets of assump- 
tions have been used in the published literature to specify the 
inner and the outer boundary of a HZ. Instead of considering ad- 
vantages and disadvantages of various definitions we use empir- 
ical boundary limits deduced from the observations of our own 
solar system. We know that the Earth is within the HZ of our 
Sun, while Venus is outside the inner boundary, and Mars is out- 
side the outer boundary. 

At some point in early stages of planetary evolution all three 
considered planets (Venus, Earth and Mars) were stripped of 
their original hydrogen rich atmospheres. The secondary atmo- 
spheres were built up by releasing mostly water vapor and car- 
bon dioxide from planetary interiors. Due to the differences in 
planetary effective surface temperature, gravity, kinematics and 
geological composition, the atmospheres of the three considered 
planets, evolved into the current diverse environments. 

Although Venus, Earth and Mars have very different surface 
albedos today (0.6 for Venus, 0.3 for Earth and 0.15 for Mars) 
we can envision that at some point in the past, after the orig- 
inal atmosphere was lost and before the secondary atmosphere 
was formed by interior out-gassing, all three planets had approx- 
imately the same surface albedo given basically by surface re- 
flectivity of the rocky surface material. We assume that the sur- 
face albedo at that time was around 0.1, somewhere in between 
the current albedo of Mars (a = 0.15) and Mercury (a = 0.05). 
To obtain effective planetary surface temperatures at the time af- 
ter the planetary primitive atmosphere was lost and before the 
secondary atmosphere was formed, we use a simple expression 
for an effective planetary temperature, T e //, implied by a top-of- 
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Fig. 1. Top panel: Effective temperatures of Venus (V), Earth (E) 
and Mars (M) as a function of planetary albedo. Bottom panel: 
Effective temperatures of "naked" V, E, M, planets c and d af- 
ter the loss of their primitive atmosphere and before acquiring 
their secondary atmosphere, and empirical boundaries of a HZ 
(horizontal red and blue lines). 



atmosphere radiative balance, 



T r S (l a), 1/4 

leff ~ i — J 



4cr 



(1) 



where, S is the incident solar flux at considered location and 
time, a is planetary albedo, and cr the Stefan-Boltzman con- 
stant, and we assume the planet to be a black body radiator 
(ei = 1) and a gray absorber (fy = 1 - a). We also assume 
the solar luminosity to be reduced by 26% (Harwit 2006 - for 
about 4.5 billion years ago) and we calculate the T e fj as a func- 
tion of the planetary albedo. Considering the surface albedo to 
be somewhere in between 0.05 and 0.15 (current surface albe- 
dos of Mercury and Mars), we detect in Fig. 1: top panel, that 
the effective temperatures were around T e ff(Mars) = 200 A, 
T e ff(Earth) = 25QK and T e ff(Venus) = 290K. From empiri- 
cal observations we know that of three planets only Earth has 
made it into the HZ. Although it is difficult to determine the in- 
ner and the outer boundaries of the HZ, we take the half points 
between Earth and the other planet's (Venus for the inner and 
Mars for the outer boundary) effective temperature as the bound- 
ary of Earth's HZ (Fig. 1: bottom panel). We consider all three 
planets to have at the time a surface albedo of 0.1, and we assign 
half of the HZ temperature width as the uncertainty of our esti- 
mate. Thus we obtain the lower temperature limit of the HZ at 
the time before formation of the secondary planetary atmosphere 
to be Tmin - 225 + 14 K and Tmax = 270 + 11 K for the upper 
limit. Translating the HZ boundary within our solar system into 
the distance from the Sun, we determine Rmin - 0.82 + 0.09 AU 
andR MAX = 1.19 ± 0.09 AU. 



3. G581 c and d planets 

Now we consider the question of whether it is possible for plan- 
ets c and d to develop into a stage in which they can keep liq- 
uid water on their surface. We assume that the primitive original 
atmospheres were swept away by stellar flares during an early 
evolution of G581 and we further assume that its luminosity did 
not change significantly after its arrival on the main sequence. 
After the original primitive atmosphere of the planets was lost 
and before a possible secondary atmosphere could be built there 
was a time when the planets were without any atmosphere. The 
albedo of these planets was solely determined by reflectivity of 
its surface rocky layer. We, therefore, assume the planetary sur- 
face albedo of 0. 1 and calculate the effective surface temperature 
of the planets c and d during the time after they lost their prim- 
itive atmospheres and before the formation of the secondary at- 
mospheres by out-gassing from the planets' interiors. 

The results (Fig. 1: bottom panel) indicate that both planets 
are outside the range of temperature that would allow a HZ with 
a possibility of liquid water on the surface. It is apparent that 
the temperature of the planet c is too high (higher than tempera- 
ture of Venus before the runaway greenhouse effect) to develop 
a habitable environment and the temperature of planet d is too 
cold (colder than the temperature of Mars). 

We suggest that under conditions similar to our solar system, 
planet c would proceed on a path similar to Venus and would be- 
come too hot to allow a liquid water to exist on its surface. On 
the other hand planet d would follow a path similar to Mars and 
would stay too cold. We note that during the considered situa- 
tion, namely before the formation of the secondary atmosphere, 
planet c was warmer than Venus by at least 35K and planet d 
was colder than Mars by at least 15 A". 

We conclude that neither of these planets could have devel- 
oped conditions allowing a liquid water to exist. This conclu- 
sion is subject to the assumptions that the processes involved in 
the evolution of the secondary atmospheres proceeded on these 
planets in a fashion similar to the processes that took place in 
our solar system within the Venus, Earth and Mars environment. 
Next we consider conditions that might have been sufficiently 
different from conditions in the solar system to allow a colder 
environment for planet c and a warmer environment on planet d. 



4. Planetary climate 

The secondary atmosphere, produced by out-gassing from the 
planetary interior, is composed mainly of water vapor and carbon 
dioxide. In the case of a suitable Earth-like environment, most of 
the water vapor condenses to form oceans and ice caps with only 
a minor part staying in the atmosphere. Carbon dioxide becomes 
bound to rocks and absorbed by oceans and the biosphere, again 
leaving only a minor portion in the atmosphere. 

The formation of a secondary atmosphere has two major 
means to affect the surface temperature. The first one is through 
the "greenhouse effect" which essentially reduced the flux of 
thermal infrared radiation to space. The second mode of sur- 
face temperature modification is through the planetary albedo 
changes connected predominantly to the cloud formation. In the 
case of Earth, the planetary albedo is decreased (compared to 
the bare rocky surface) by the presence of oceans, and increased 
by the clouds in such a way that the average planetary albedo 
changes from its estimated original "rocky surface" value of 0. 1 
to the current Earth planetary albedo of about 0.30. The albedo 
change by itself would decrease the effective radiative temper- 
ature of the Earth by about 16K, from the value of 271 K (for 
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Fig. 2. Effective temperatures of "naked" Venus (V), Earth (E) 
and Mars (M), boundary of HZ (horizontal green and red lines), 
and effective temperatures of planets c and d as a function of 
their planetary albedo. 



the current value of solar flux and albedo of 0.10) to 255K (the 
current value of solar flux and albedo 0.30). 

However, the water vapor and carbon dioxide greenhouse ef- 
fect contributes to the warming of about 33K and increase the 
average of surface temperature to about 288/T. About a half of 
greenhouse warming is utilized to compensate for increase of 
albedo (from 0.1 to 0.3 and connected temperature decrease of 
16K) and the other half produces an increase of the surface tem- 
perature, by about 11 K, from the original 21 IK to 288ZT. 

If the temperature of a planet is too low, the water vapor and 
carbon dioxide will condense into ice and will be removed from 
the atmosphere and deposited on the surface. The result is an in- 
creased planetary albedo and an atmosphere lacking any green- 
house effect (Mars-like environment). In the opposite case of too 
high a planetary temperature, all water vapor and carbon dioxide 
stay in the atmosphere. The planetary albedo will significantly 
increase due to a high reflectivity of water or CO2 clouds, how- 
ever, the strong greenhouse effect will dominate and compen- 
sate for an increased albedo and cause an additional significant 
warming of the planet (Venus-like environment). The situation 
of planets c and d is depicted in Fig. 2. To bring planet c into a 
HZ, it needs to increase its albedo from the bare planet value of 
0. 10 to about 0.69 to 0.70. This can hardly be done without cov- 
ering the whole planet by clouds, which in turn implies a dense 
H2O vapor and/or CO2 atmosphere. Such an atmosphere would, 
however, create a large greenhouse warming that would push the 
planet back to the region of too high temperature for liquid water 
to exist on its surface. Thus it seems that under the considered 
scenario the planet c has a little chance of developing conditions 
for liquid water to exist on its surface. 

The bare planet d is too cold to hold liquid water (T e ff = 
183AT). We can assume that a medium strength (comparable 
to Earth) atmospheric greenhouse effect could provide about 
35K warming. The accompanying cloudiness would increase the 
planetary albedo to about 0.30, which would produce a cooling 
of about II K. Thus the result of Earth-like atmosphere would 
produce a net surface warming of about 24K, which is not suffi- 
cient to bring the planet into the range of habitability. However, 
a considerably denser atmosphere that would lead to a planet 
completely covered by clouds, producing a planetary albedo of 
about 0.6, and warming of around 100K (about three times the 
greenhouse effect of the Earth's atmosphere) would be sufficient 
to bring the planet d into a HZ. No objectionable mechanism 
is apparent that would prohibit the formation of such an atmo- 
sphere. 



We conclude that assuming a situation similar to our solar 
system, it is very unlikely that the planet c could be within a HZ. 
Whereas, the planet d with an atmosphere producing three times 
the greenhouse effect on Earth, could be within the HZ. 

5. Atmospheric circulation on a gravitationally 
locked planet 

Climate and atmospheric circulation on gravitationally locked 
planets is not currently well understood. For the case of G581 
planets, based on kinematics considerations, von Bloh et al. 
(2007a) invoked the notion that planets inside the HZ of M stars 
are tidally locked. To better understand and justify this assump- 
tion we tested this hypothesis in our solar system planets. 

We define the dimensionless parameter, <f>, as the ratio be- 
tween the planet orbital period (in terrestrial days) to the length 
of the planetary day (in terrestrial hours). Obviously, for tidally 
locked planets, (p is equal to 0. In Table 1, we present the values 
of <p for the eight recognized planets noting that at the distance 
of Mercury from the Sun, d = 0.39 AU, (f> value is small sug- 
gesting that Mercury is close to being a tidally locked planet. 
Consequently, at the estimated distances of planets c and d from 
the central star G581, of 0.073 AU and 0.25 AU, respectively, 
the assumption of gravitationally locked planets seems reason- 
able when compared with our solar system. We also included 
in Table 1 for G581 and solar system planets the gravitational 
strength values demonstrating that planets c and d are consider- 
ably more tidally locked than any planet in our solar system. 

Now we present basic qualitative ideas and we emphasize the 
need to modify existing atmospheric general circulation models 
to investigate the problem in detail. Since only a half of a planet 
is exposed to incident solar radiation, the radiative energy bal- 
ance equation needs to be modified to reflect this situation. We 
have, under similar assumptions as in (1), 



S (l - a)nr 2 = 27rr 2 o-(T 4 + TV), 



(2) 



where, S is an incident radiative flux, a is an albedo of an il- 
luminated part of a planet, r is planets radius and T and T a are 
the effective temperatures of an illuminated and shadow side, 
respectively. If there is no energy transport between the illumi- 
nated and shadow side of a planet, or if the transport is negligible 
so that T a « T can be neglected we have, 



l eff 



Soil 



' a )-il/4 



la 



-]' 



(3) 



and the effective temperature on the illuminated side increases 
by about 19% compared to the case of planet that is not grav- 
itationally locked-in. At the same time the temperature on the 
shadow side decreases and without atmospheric circulation and 
possible heat conduction through the planetary interior it would 
be close to the cosmic background temperature. However, the 
differential heating of the illuminated side will produce atmo- 
spheric circulation, which will initiate a heat transport between 
the illuminated and not illuminated hemisphere. The maximum 
rate of heating will occur at the center of the illuminated side 
of the planet. Here the air rises up creating a low pressure region 
near the surface of the central region of the illuminated side (Fig. 
3). We acknowledge that serious climate stability considerations 
are problematic under this scenario (Porto de Mello, del Peloso 
& Ghezzi 2006). 

Due to the axial symmetry of the illumination the surface air 
from all directions will move towards the low pressure at the cen- 
ter of the illuminated side while the upper air will diverge from 
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Table 1. Solar system values of <p and the value of the pull or gravitational strength (kg 2 /m 2 ) relative to Earth. 
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Fig. 3. Conceptual illustration of atmospheric circulation on a 
gravitationally locked-in planet. 

the center into all directions towards the edges of the shadow 
side. The upper air will experience radiative cooling as it moves 
from the center towards the edges and, when sufficiently cool, 
will start descending towards the surface. This thermally driven 
circulation is an analogy to the Hadley circulation cell of the 
Earth atmosphere. The differences are, however, quite signifi- 
cant (e.g, the absence of a Coriollis force and an axial symmetry 
of circulation on a gravitationally locked-in planet). 

In this way, two distinct regions of the interaction between 
the air of the illuminated and shadow sides can develop. Still 
partially warm upper air moving away from the central region of 
an illuminated side can "overshoot" and enter the shadow part 
of the planet. The low level cold surface air from the shadowed 
side can flow into the lower pressure regions towards the center 
part of the illuminated side. How much of these possible air ex- 
changes really takes place will depend on the parameters of the 
locked-in planet atmosphere that enters into the Navier-Stokes 
equations of motion. 

Qualitatively the atmospheric circulation will try to decrease 
the temperature gradient between illuminated and shadow parts 
of the planet. To assess how much heat could be transported if 
planets c and d are gravitationally locked-in will require detailed 
knowledge of the atmospheric properties of the planets. A sim- 
ple example based on an energy balance indicates that to keep 
the shadowed part at the temperature of 2007^ and the illumi- 
nated part at 3QQK would require about 20% ((T () /T) 4 ~ 0.20) 
of the incident energy to be transported towards the shadow side 
of the planet, while 25QK and 300K temperatures would require 
almost 50% ((T /T) 4 ~ 0.48) of incident energy to be trans- 
ported to the shadowed side. Thus it would be extremely dif- 
ficult to maintain on the planet d the shadow side temperature 
at the level that would prevent water vapor and carbon dioxide 
from being removed from the atmosphere. On the other hand 
planet c, in the locked-in situation may have a chance to develop 
a habitable region. To keep the shadowed site of the planet at the 
temperature of about 250K would require about 26% of the heat 
to be transported to the shadowed side, while about 15% trans- 
port would be sufficient to maintain the average temperature at 
22QK. Considering a possible greenhouse effect of about double 
that found currently in Earth's atmosphere, we can reach a sur- 



face temperature around 2867^, just right to almost duplicate our 
earthly environment. Of course, the problem would be the over- 
heated illuminated side of the planet at about 35QK and a lack 
of sunshine on the shadow side. Perhaps tropical forests on the 
hot side can provide a sufficient amount of oxygen for the whole 
planet and near the illuminated/shadowed boundary a habitable 
region may occur. 

6. Conclusion 

We use an empirically defined boundary of a HZ in analogy with 
our solar system and the fact that only the Earth is within the HZ, 
while Mars and Venus are not. By assuming the formation of 
the secondary atmosphere similar to Earth, we find that neither 
planet, c nor d, would be within a HZ. However, for planet d 
there is a possibility to form a HZ provided that a greenhouse 
effect of about 100K could be developed. Conversely, planet c 
seems to be too hot and we are unable to foresee any situation 
under which a region of habitability could be developed, under 
Earth-like conditions. 

If the planets are gravitationally locked-in, as we have shown 
to be a reasonable assumption, there seems to be a possibility 
especially for planet c to develop atmospheric circulation that 
would transport heat from the illuminated to the shadow part of 
the planet and form a definite region on the planet where liq- 
uid water could exist. The future assessment of this possibility 
will require the development of a simplified Budyko-style cli- 
mate model (Budyko 1969, Chylek and Coakley 1975), as well 
as modification of atmospheric general circulation models and 
additional information concerning the planets and their atmo- 
spheric composition. 

We found that a gravitationally locked-in planet can have liq- 
uid water on its surface even if it is outside a HZ, defined tra- 
ditionally by the parameters and distance to the sustaining star 
(Udry et al. 2007). Since half of the planet is permanently in 
"sunshine" an atmospheric circulation can develop between the 
illuminated and shadow parts and create a HZ with a possibil- 
ity of surface liquid water in a planet that, otherwise, would be 
unsuitable. Thus the auxiliary definition of a HZ should be ex- 
tended to include the rotational and atmospheric characteristics 
of planets; a concept alluded to by von Bloh et al. (2007a, b). 
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